Wave packet evolution in non-Hermitian quantum systems 
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The quantum evolution of the Wigner function for Gaussian wave packets generated by a non- 
Hermitian Hamiltonian is investigated. In the semiclassical limit h — »> this yields the non-Hermitian 
analog of the Ehrenfest theorem for the dynamics of observable expectation values. The lack of 
Hermiticity reveals the importance of the complex structure on the classical phase space: The 
resulting equations of motion are coupled to an equation of motion for the phase space metric — a 
phenomenon having no analog in Hermit ian theories. 



Effective non-Hermitian Hamiltonians have long been 
used for the description of a wide range of open quan- 
tum systems [1 j. Their applications range from chemical 
reactions to ultra cold atoms and laser physics. Com- 
plex potentials for matter waves can be tailored exper- 
imentally using standing light waves [2], and the com- 
plex Schrodinger equation appears in optics using ma- 
terials with complex refractive index. The latter anal- 
ogy was used in the recently reported first experimental 
realizations of PT-symmetric non-Hermitian Hamiltoni- 
ans [31 H], which boosted the interest in the field further. 
While much attention has been paid to the theoretical 
study of example systems [5j [6], the generic dynamical 
features remain hitherto mostly unexplored. The present 
paper aims to fill this gap by investigating wave packet 
dynamics for general non-Hermitian systems. 

The time evolution of wave packets is a powerful tool 
for the understanding of both dynamical and stationary 
properties of Hermitian quantum systems pj [8] . Further- 
more, it is a convenient way to investigate the semiclas- 
sical limit and thus forms the basis of many semiclassical 
methods. In what follows we shall generalize the funda- 
mental ideas of wave packet dynamics to non-Hermitian 
systems. From this we derive classical equations of mo- 
tion in the spirit of a generalized Ehrenfest theorem. This 
will pave the way for the development of a semiclassical 
framework for non-Hermitian quantum systems, or more 
generally, absorbing wave equations. 

The semiclassical properties of non-Hermitian dynam- 
ics have recently been approached from various direc- 
tions. Examples include the study of ray dynamics of 
absorbing wave equations for weak non-Hermiticities [9]; 
the mean-field approximation for a non-Hermitian many- 
particle system [6] ; the quantum classical correspondence 
for open quantum maps in the chaotic regime [TP] : and 
complex extensions of the quantum probability distribu- 
tion [11 . In [12] a coherent state approximation has 
been applied to non-Hermitian systems, and a general- 
ized canonical structure involving a metric gradient flow 
has been identified. Here we go beyond this study by in- 
vestigating general Gaussian states that are allowed to 
change their shapes during time evolution, which can 
be interpreted as a time-dependent metric on the cor- 



responding classical phase space. To derive the classical 
evolution equations in the spirit of the Ehrenfest theo- 
rem, we study the quantum evolution equation for the 
Wigner function and take the semiclassical limit. This 
results in a new type of classical phase space dynam- 
ics in which the evolution equations for the phase space 
coordinates depend on the local metric, and vice versa. 
The complex structure of the classical phase space, which 
is rarely considered in Hermitian systems, thus becomes 
highly relevant in the presence of non-Hermiticity. The 
main dynamical effect of the anti-Hermitian part of the 
quantum Hamiltonian in this semiclassical approxima- 
tion is a damping of the motion. This strengthens the 
often speculated connection to classical dissipation where 
non-Hermitian Hamiltonians are a recurrent theme in the 
search for quantum counterparts (see e.g., [13 and refer- 
ences therein). 

Gaussian coherent states. For a general n-dimensional 
quantum system we study the evolution of initial wave 
packets of the form 

_ (det ImP) 1 / 4 u P . {q _ Q)+ i {q _ Q) . B{q _ Q)] 

where P, Q G R n , and B is a complex symmetric ma- 
trix with positive definite imaginary part so that ip is 
localized around q — Q and normalized to unity, he 
Gaussian states ([!]) with fixed B form a submanifold 
of Hilberts space wich has a natural complex structure 
associated with B. In the semiclassical limit this sub- 
manifold can be identified with the classical phase space, 
which thus inherits not only the symplectic but also a 
metric structure. These relations become most trans- 
parent using a phase space representation [14]. In the 
following we will focus on the evolution of the Wigner 
function W of ^. The Wigner function alllows a direct 
computation of expectation values via phase space inte- 
grals (^i^)/^,^) = (A) w := JWAdz/ JWdz, where 
A is the Weyl quantization of A(z) and z = (p, q) are 
canonical phase space coordinates. The Wigner function 
of a Gaussian state ([I]) is Gaussian: 

W(z) = (nh)- n e-^ z - z ^ z - z \ (2) 
Here Z = (P, Q) e R n x R n and the matrix G is related 
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The matrix G is nondegenerate, positive, and symmet- 
ric, and thus acts as a metric on phase space. It is also 
symplectic, i.e. it satisfies GQG = ^, where 



-I 

1 



(3) 



is the symplectic, or canonical, structure on phase space. 
With such a metric we can associate a compatible com- 
plex structure J with J 2 = —I and VtJ = G [15 . The 
structure defined by £7, G, and J turns the phase space 
into a Kahler manifold. 

The Wigner function ([2| is localized of order \fh 
around the maximum Z = (P, Q). In the semiclassi- 
cal limit H —> this collapses to a phase space point Z. 
Hence the expectation value of an observable A satisfies 



(A) w = A(Z) + 0(h) 



(4) 



for A(z) smooth. The metric G describes the shape and 
orientation of W in phase space, and therefore determines 
the variance of observables: 



h 



VA(Z)-G- 1 VA(Z) + 0(h 2 ) 



(5) 



In the Hermitian case an initially Gaussian state stays 
approximately Gaussian during the time evolution up to 
the Ehrenfest time [7]. The center moves according to 
the classical canonical equations of motion Z = Q\7H, 
and the evolution of the metric is governed by the lin- 
earized Hamiltonian flow around the classical trajectory. 
In the framework of the time-dependent variational prin- 
ciple [16 it has been shown that this dynamics can also 
be described by Hamiltonian equations of motion. Al- 
though it plays a central role in semiclassical methods 
involving families of coherent states [I6j [17], the met- 
ric is often little investigated, as it does not enter the 
dynamical equations for the phase space variables. We 
will see shortly that this is fundamentally changed in the 
presence of non-Hermiticity. 

Non- Hermitian Wigner-von Neumann equation. De- 
composing the Hamiltonian in its Hermitian and anti- 
Hermitian part H — if, where we assume H and T to 
be given as the Weyl quantizations of sufficiently well- 
behaved classical observables H(z) and r(z), the evolu- 
tion equation for a density operator W follows from the 
Schrodinger equation as 



ihd t W=[H,W]-i[f,W} + , 



(6) 



where [•,•]+ denotes the anti commutator. Thus, the evo- 
lution equation of a general Wigner function is given by 

ihdtW = (HftW - WftH) - i(IW + W$T) , (7) 



where (A$B)(z) denotes the Weyl product [14] for two 
phase space functions A{z) and B(z): 

A$B = Ae^-- a *-B ~jr^ (| ) k A(^ z ■ ^ z ) k B . 

k=0 ' ^ ' 

Here the arrows over the differential operators indicate 
whether they act on the function to the right or to the 
left. We will now evaluate the leading order terms in 
H of 0. The Hermitian part of the evolution equation 
Q is the well known Moyal bracket with an asymptotic 
expansion in odd powers of H whose leading term gives 
the Poisson bracket 



HftW - W$H = ihVH ■ nvw + 0(h 3 



(8) 



The anti-Hermitian part has an expansion in even powers 
of h, with the first two terms given by 

IW + W$T = 2TW - —A r W + 0(ft 4 ), (9) 

where we introduced a second order differential operator 
defined by T as A r W := r(K7 z ■ z fw. Denoting the 
matrix of second derivatives of T(z) at z by T"(z) the 
operator Ar can be written in the form Ar = V • T^V, 
with := VL l Y"(z)VL. It follows that A F is Hermitian. 

Furthermore, \iT"{z) is symplectic, then Tq = and 
Ar is the Laplace-Beltrami operator defined by T". 

Summarizing, in leading order of h the dynamical equa- 
tion for the Wigner function reads 



- A F - hVH • OV + 2Y )W . (10) 



hd t W = 



For vanishing T we recover the classical Liouville equa- 
tion for the transport of phase space densities. For non- 
vanishing and positive T, on the other hand, the T term 
defines a diffusion equation. 

The higher order terms are of order h 3 \d 3 W\(\d 3 H\ + 
|9 3 r|), i.e. they are small if the derivatives of H(z), 
T(z) and W(z) do not grow too fast as h —> 0. If the 
derivatives of W are bounded, the first term on the right 
side of ( p~Q| ) is of lower order than the other terms. In 
this case the solution W(t, z) is obtained by transporting 
the initial W(z) along the Hamiltonian flow generated by 
H, multiplied by a damping factor, which is determined 
by the integral of T along the Hamiltonian trajectories 
of H. This behavior is well known from damped wave 
equations. For a Gaussian initial state (J2|, on the other 
hand, the term H 2 AyW = 0(h) is of the same order as 
the Hamiltonian term in ( [To] ), and the dynamics differ 
drastically from the Hermitian case. 

Gaussian evolution. In what follows we investigate the 



solution of (10) for an initial Gaussian Wigner function 
([2]). Inserting a Gaussian ansatz for the time evolved 
Wigner function 



W(t,z) = 



^_ e -i5z.G W Sz with §Z :=z- Z(t) 



3 




into ( |10 ) yields 

a 



H-+2Z- G5z - 5z • G5z 
a 



W(z) 



Sz ■ GT'ASz - 2VH ■ QGSz 



tr (r^G)-2r 



(11) 



W(z) 



Following the well established method of Heller and Hepp 
[7], we expand T(z) and H(z) up to second order around 
z = Z: T(z) « T(Z) + Vr(Z) • Sz + \5z ■ T"(Z)5z and 
Vff (z) « Vff(Z) + H"{Z)5z. Since is localized 

around z = Z with a width proportional to the re- 
mainder terms are of order ft 3 / 2 . Separating different 
powers of Sz = z — Z in ( 11 ) then yields the following 



three equations of motion for G(t) and a(t): 

z = n\/H(z) - G~ 1 vr(z) (12) 

G = H"(Z)QG - GQH"(Z) + T"(Z) - GT&(Z)G (13) 



a 
a 



1 



r(z)--tr[rS(z)G] . 



(14) 



To obtain ([13]) the symmetry enforcing convention G = 
(G £ + G)/2 was applied. As W depends only on the sym- 
metric part of G any anti-symmetric part is unobservable. 

The time evolution of expectation values and vari- 
ances of arbitrary observables in Gaussian coherent states 
for small h is determined by Z(t) and G(t) according 
to Q and ([5]). Equations (fl2|)-([T4|) can be interpreted 
as the hitherto unidentified semiclassical limit of non- 
Hermit ian quantum dynamics. This result goes beyond 
previous studies of the non-Hermitian Ehrenfest theo- 
rem p~2j [13j [18] for two reasons. First, previous stud- 
ies usually focussed on unnormalised expectation values, 
which prevented the identification of a classical structure, 
and second, disregarded the role of the metric, related to 
the widths of the quantum wave packet. The dynamics 



(12) emerging as the classical limit is no longer Hamil- 



tonian, but has a Hamiltonian part and a gradient part, 
determined by the Hermitian and anti-Hermitian parts 
of H — if, respectively. The main dynamical effect of the 
anti-Hermitian part is to drive the motion towards the 
minima of T. In addition, this gradient part is coupled 



to an evolution equation (13) for the metric G which in 



turn depends on (12). In this context it is important 



to note that (13) preserves the symplectic nature of G 



and hence describes an evolution of the complex struc- 
ture on phase space. Further, the anti-Hermitian part 
leads to a change of the overall probability according to 



equation (14), which can be interpreted as absorption or 



amplification. The first term gives the contribution from 
the center and the second term captures the influence of 
the width of the Wigner function. Note, however, that af- 
ter renormalization the non-Hermitian Schrodinger equa- 
tion is equivalent to norm-conserving nonlinear models 
for quantum dissipation [19] • 






FIG. 1. Time evolution of the exact Wigner function (left 
column) and the semiclassical approximation (right column) 
for an initial state at (p, q) = (5, 0) at different times (t = 
1, 2.5, 4) for the anharmonic oscillator. The white line shows 
the motion of the center. The left panel on the bottom 
shows the norm of the exact quantum state (black dashed 
line) and the semiclassical approximation (blue line), and the 
right panel shows the largest eigenvalue of G(t) (blue line) in 
comparison with the Hermitian case 7 = (pink line). 



The quadratic approximation around z = Z(t) to H(z) 
and T(z) is expected to remain accurate so long as W(t, z) 
stays strongly localized around z = Z(t). Since for a 
symplectic G we have ||G _1 || = ||G||, a suitable criterion 
for this is 



(15) 



The wave packet becomes delocalized at the Ehrenfest 
time Te defined by S||G(Te;)|| = 1 and the semiclassi- 
cal approximation based on the central trajectory Z(t) 
breaks down. The nonlinear term in the equation for 
G(t) that is induced by T can have a stabilizing effect 
on the long-time evolution of G(t). Therefore, the non- 
Hermitian part can increase the Ehrenfest time, i.e. the 
time scale for which the semiclassical approximation is 
valid as compared to the Hermitian case. 

Examples. To illustrate our results we consider two 
examples. The first example is a non-Hermitian an- 
harmonic oscillator with H = ^(p 2 + q 2 ) + f<? 4 and 
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FIG. 2. Quantum evolution (black dashed line) versus semi- 
classical approximation (blue line) of a PT-symmetric waveg- 
uide for an initial state at (p, q) = (0, 2). Shown are the phase 
space evolution (left), and the evolution of the norm (right). 



tern. 

Conclusion. The results presented here for the evolu- 
tion of a Gaussian coherent state generated by a non- 
Hermitian Hamiltonian provide the basis for a more pro- 
found understanding of non-Hermitian time evolution, a 
topic of considerable interest in a wide range of subjects. 
In particular, the application to realistic examples of typ- 
ically non local non-Hermitian Hamiltonians appearing in 
resonance physics is an interesting task for future studies. 

The authors thank D. C. Brody and H. J. Korsch for 
useful comments. EMG is supported by an Imperial Col- 
lege Junior Research Fellowship. 



f = ^(p 2 + q 2 ) respectively, with uj = 1, 7 = 0.2, and 
j3 = 0.5. This can be interpreted as a quantum analogue 
of a damped anharmonic oscillator (see, e.g., [13 and 
references therein). The simple structure of this model 
makes it an ideal testing ground for the semiclassical 
approximation proposed here. We set ft = 1, which is 
equivalent to a rescaling upon which j3 plays the role of 
an effective ft. For j3 = the semiclassical approxima- 
tion becomes exact. While the Hermitian part tries to 
propagate a state along closed curves of constant ener- 
gies around the origin, the anti-Hermitian part drives it 
towards the origin, thus acting as a damping. Figure [I] 
shows the exact numerical propagation and the semiclas- 
sical approximation for an initial Wigner function with 
Go = which are in good agreement. Also the total 
mass of the exact state, a measure for the absorption, 
due to the anti-Hermitian part is well described by the 
semiclassical approximation a(t), as illustrated in the left 
panel on the bottom. The right panel shows the time- 
dependence of the larger eigenvalue of G as a measure for 
||G(t)|| in comparison with a Hermitian case 7 = 0. The 



result indicates that the Ehrenfest time, see (15), is in- 



creased in the non-Hermitian case, i.e. the semiclassical 
approximation is accurate over a longer time scale than 
in a comparable Hermitian case. 

Second, we consider a simple model system for a PT- 
symmetric optical waveguide [3] : A single waveguide with 
harmonic confinement described by H = ^(p 2 + g 2 ), and 
an anti-Hermitian part T = 5tanh(0.2g) that models ab- 
sorption on one side and equally strong amplification on 
the other side with a smooth transition in between. Al- 
though the Hamiltonian is complex, due to its special 
symmetry, it has real eigenvalues which can lead to a 
pseudo-closed behavior. This phenomenon is captured by 
our classical approximation. Figure [2] shows an example 
of the full quantum evolution and its classical counter- 
part. Both the phase space evolution and the dynamics 
of the norm are well approximated by the classical de- 
scription. In particular, despite the anti-Hermitian part 
in the Hamiltonian, no sink of the dynamics is observed. 
Note that there is a stable fixed point at (p, q) = (0, 1), 
corresponding to the ground state of the quantum sys- 
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